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Abstract 

Angular- momentum-projected energy surface calculations for A ~ 110 nuclei indi- 
cate three distinct energy minima occurring at different angular-momenta. These 
correspond to normal, super-, and hyper-deformed shapes coexisting in one nu- 
cleus. 110 Pd is studied in detail, with a quantitative prediction on super- and hyper- 
deformed spectra by the Projected Shell Model calculation. It is suggested that 
several other neighboring nuclei in the A-110 mass region, with the neutron number 
around 64, also exhibit clear super- and hyper-deformation minima. 
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1 Introduction 

The usual assumption of nuclear structure models is that a nucleus can be de- 
scribed as a set of interacting particles moving in an average potential. The nu- 
cleonic motion in a potential well often occurs in a highly coherent way, giving 
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rise to collective phenomena such as rotational spectra and large quadrupole 
moments. These led to the idea of introduction of the concept of nuclear defor- 
mation [Tj. Since then, the search for exotic nuclear shapes has continuously 
been one of the research frontiers in nuclear structure physics. One remark- 
able example is the discovery of super-deformed (SD) rotational bands [2]. 
Moreover, microscopic calculations [3,4,5f6] have predicted that some nuclei 
may attain even greater quadrupole deformations corresponding to a hyper- 
deformed (HD) shape, and that the nuclei with these deformations are possi- 
bly stable against nuclear fission. Experimentally, the search for HD states has 
been continued for years (see, for example, Ref. [7] and conference proceedings 
of the 2004 Zakopane School [8]). 



The observation of rotational bands in 108 Cd has stimulated a new wave of 
study along this line. In Ref. [9], Clark et al. reported a rotational band with 
the data indicating that this nucleus has the most deformed shape identified 
to date. Using the Projected Shell Model (PSM) [10] , some of us analyzed 
the structure of the band. It was found [IT] that the wavefunctions of the 
observed states notably contain a component of the proton i 13 / 2 orbital from 
the N = 6 harmonic oscillator shell. It is well-known that the 7ri 13 / 2 orbitals 
usually appear near the proton Fermi level in normal- deformed (ND) nuclei 
of the actinide region, and near those in SD nuclei of the rare-earth region. In 
a Z = 48 nucleus of the A-110 mass region, occupation of the 7u 13 / 2 orbitals 
is possible only if the nucleus has a greater deformation. 



Inspired by these advances, we have employed the PSM to carry out a system- 
atical study on the neighboring nuclei of 108 Cd. Angular-momentum-projected 
energy surface calculations have explored distinct minima in some of these nu- 
clei, separately lying at different deformations and angular momenta. Based 
on the content of the wavefunctions, in which nucleons from several major 
shells contribute to the collectivity near the Fermi surface, we can understand 
the occurrence of the minima as a correlation effect involving multiple major 
shells. According to our analysis, 108 Cd has two energy minima. The observed 
band in 108 Cd [H] (with neutron number 60) belongs to the second minimum. 
As we shall show in the present paper, three energy minima, which correspond 
to ND, SD, and HD shapes coexisting in one nucleus, can be formed in nuclei 
with neutron number around 64. 



The paper is organized as follows. In Section 2, we outline the theory and 
give the calculation conditions. In Section 3, we take 110 Pd as an example to 
present the main results of the present study. Further quantitative predictions 
are given in Section 4. Finally, the paper is summarized in Section 5. 
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2 Theory of angular momentum projection 

Using a conventional shell model to study deformed, heavy nuclei is desirable. 
However, it is unfeasible in practice because of the unavoidable problem related 
to huge basis dimensions. Therefore, the study of deformed, heavy nuclei has 
relied mainly on the mean field approximation. In a mean field theory, single 
particle motion is treated properly which gives rise to a correct description of 
shell structure. However, correlations beyond the mean field level are usually 
neglected. If the energy gain due to the correlations would be a constant with 
no particle-number and angular-momentum dependence, then consideration of 
correlations may not be important. However, it is often not the case in reality. 
As we shall see in the 110 Pd example, when the correlation effect is added to 
that of shell structure, several energy minima clearly occur in the projected 
energy curves. 

The Projected Shell Model [TO] is a shell model built from a deformed basis. To 
be concrete, the PSM constructs its shell-model space by using the deformed 
Nilsson single-particle states with a deformation 62- Pairing correlations are 
incorporated into the Nilsson states by the BCS calculations. The consequence 
of the Nilsson-BCS calculations defines a set of quasiparticle (qp) states as- 
sociated with the qp vacuum |0(e 2 )) = |0) in the intrinsic frame. One then 
builds shell model bases by writing multi-qp states on top of |0): 



where ot* is the qp creation operator and the index n (p) denotes neutron 
(proton) Nilsson quantum numbers of the orbitals. In Eq. ([]]), when the higher 
order multi-qp states denoted by "• • ■" are all included, one recovers the com- 
plete shell model configuration. In practice, these high order multi-qp states 
are not necessarily needed if one is interested in the near-Yrast property only. 
We have checked that for producing energy minima reported in the present pa- 
per, nearly the same results are obtained when the basis size in ([T]) is changed 
from a few to a few hundreds. Angular momentum projection transforms the 
wavefunction from the intrinsic frame to the laboratory frame. The total wave- 
function is expressed as 



where P 1 is the angular-momentum-projection operator [12] and k labels the 
basis states. Finally a two-body shell model Hamiltonian is diagonalized in the 
projected space, and the diagonalization determines fj. in This relatively- 
simple model has been extensively tested in the past decade. For the applica- 
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tions related to the present work, the PSM has been successful in describing 
the SD states in different mass regions [T3|14fl5lfT6] . 

In searching for energy minima, we calculate the angular-momentum-projected 
energy 



E 1 





H 









(3) 



as a function of static quadrupole moment defined by 



W> = e °» </020| JO) • 



where the effective charges are e c g = 1.5e for protons and 0.5e for neutrons. 
The energy E 1 and the variable (Q) are both evaluated by using the same 
wavefunction (T5]). 

Equation (3) represents a study of energy minimization after angular mo- 

\E ri ) can have 



mentum projection. In principle, the energies calculated with 
advantages over those with unprojected state |0) in that the projection treats 
the states fully quantum-mechanically by collecting all the energy-degenerate 
mean-field states associated with the rotational symmetry. Angular momen- 
tum projection is a superposition in the Euler space which lowers the total 
energy of the system. The energy gain is due to the correlations, which is 
particularly important in the present content for the formation of local de- 
formation minima. Although in some cases, especially when the minima are 
deep and stable, the results between projected and unprojected calculations 
are qualitatively similar, there have been many examples with qualitatively 
different conclusions (as for instance, the one shown in the seminal work of 
Hayashi, Hara, and Ring [Hj). In recent years, there have been attempts to 
incorporate angular-momentum-projection in the theories with effective in- 
teractions [T8lll9ll20|21|l22] . and notable effects due to projection have been 
found. 

The Hamiltonian appearing in ([3]) is of the quadrupole plus pairing type, with 
the inclusion of the quadrupole-pairing term [I0J 

h = h - 1 -xY: QlQ, - G M P^P -G q Y: PIP,- (5) 

t 1 /i 



This is a rotation-invariant Hamiltonian with each term in it is a scalar. In Eq. 
([5]), Hq is the spherical single-particle Hamiltonian which contains a proper 
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Table 1 

Quadrupole-quadrupole interaction strengths (in MeV). 





Xnn 


Xpp 


Xnp Xpn 


Region I 


0.0377 


0.0302 


0.0337 


Region II 


0.0355 


0.0284 


0.0318 



spin-orbit force [23]. The monopole pairing strengths Gm are taken to be 
11.5/ A for neutrons and 19.5/ A for protons, where A is the total particle num- 
ber. The quadrupole pairing strength Gq = 0.16Gm, which is the value usually 
taken in the PSM calculation. These parameters are found to be appropriate 
for a correct description of the known data. To incorporate the shell effect 
across a large single-particle space, nucleons in three major shells (N = 3, 4, 5 
for both neutrons and protons) and in four major shells (N = 4, 5, 6, 7 for 
neutrons and N = 3,4,5,6 for protons) are activated, respectively, for the 
smaller deformation region (denoted as Region I) and for the larger defor- 
mation region (denoted as Region II). The quadrupole-quadrupole interaction 
strengths \ are listed in Table I for the regions I and II, and these values are 
determined in a similar way as in Ref. [2Tj . 



3 Detailed analysis of Pd 



In Fig. 1, we show the calculated energy surfaces for 110 Pd with different 
angular momenta. The results clearly indicate three local minima appearing 
at distinct deformations (marked by red dots). The first minimum starts from 
(Q) =0.7 eb at spin 1 = 0, and as spin increases, gradually moves to a higher 
(Q) value of 1.0 (corresponding to a deformation of 62 ~ 0.25). The change of 
deformation with spin indicates a softness of the system at the low angular- 
momentum region. At about 1 = 8, the shape is stabilized and the system 
begins to show a rotor character. We note that Regan et al. [21] discussed, 
with a collection of many experimental data in this mass region, the evolution 
of the rotational behavior in terms of transition from vibrational to rotational 
motion. Our present calculation may provide a microscopic insight to the 
observation, which suggests that this rotational behavior can be understood 
as a gradual increase in deformation with increasing spin in the first energy 
minimum [25] . 

The second minimum is located at about (Q) =2.8 eb, corresponding to 62 ~ 
0.57. The minimum is fully developed only at spin I = 24 and above. The 
third one starts showing up at a rather high spin I ~ 44 at about (Q) = 
3.9 eb, corresponding to 62 ~ 0.79. Unlike the first one, the second and the 
third minimum are stable {i.e. the deformation minimum does not change 
with spin), reflecting a rigid-rotor character. It is worthwhile mentioning that 
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Fig. 1. Calculated energy surfaces in 110 Pd for different angular momenta as a 
function of static quadrupole moment (Q). (a) Region I. The energy zero is set to 
be the total energy at (Q) =£2 = and 1 = 0. (b) Region II. The absolute energy 
is arbitrarily chosen. Note that the plots do not give the relative energy between 
the two regions. 

in a sharp contrast, unprojected energies (not shown in the figure) suggest 
a qualitatively different picture: a flat bottom centering at e^, = and only 
one rather shallow minimum at 62 ~ 0.67. We believe that this difference 
comes from the treatment with projection, which incorporates the many-body 
correlations missed in the mean-field states. 

A complete description (including particle-alignment and back-bending in mo- 
ment of inertia) of the ND states usually requires a mixing of the shell with 
the next A-shell (with the high-j intruder orbitals included). Now our cal- 
culation indicates clearly that the SD and HD minima are a consequence of 
a shell-mixing with the next two and next three A-shells, respectively. This 
picture is consistent with the main conclusions of Dudek et al. |26 ; 27 ; 28j. As 
schematically shown in Fig. 2, the diagram exhibits a correspondence between 
formation of a deformed shape and an actual participation of nuclear shells in 
the collective motion. It is clear that formation of each higher level of defor- 
mation is associated with particles from an additional major shell in the shell 
mixing near the Fermi level. In our example here, the occurrence of the HD 
minimum in 110 Pd requires a collective participation of four neutron shells up 
to the N = 7, and four proton shells up to the A = 6 shell, respectively. 

For a higher A-shell, it is evident that among all the orbitals, the down- 
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Fig. 2. Schematic diagram illustrating a correlation between a formation of defor- 
mation minima and shell mixing near the Fermi level: ND <-> adding N + 1 shell, 
SD <-> adding iV + 2 shell, and HD <-> adding TV + 3 shell. 

sloping intruder orbitals with the highest j are always occupied first. It is well- 
known that these intruder orbitals are important in determining the nuclear 
high-spin behavior because of their unique properties with rotation. Thus the 
next question is what role these intruder orbitals play in the formation of the 
extreme deformation. To answer this question, we consider the expectation 
value of the occupation operator 



n 



(6) 



If the summation runs only over the highest j intruder orbitals we denote it 
as n intr (for the HD minimum, j intr is ji 5 / 2 for neutrons and Zi 3 / 2 for protons), 
otherwise as n a n when summing over all the orbitals of like particles. We 
calculate the ratio r = n intr /n a \\ for each spin state by including those orbitals 
within ±5 MeV around the Fermi level. 

In Fig. 3, the ratio r obtained at three deformations around the HD minimum 
in 110 Pd are separately shown for protons and neutrons. It can be seen that 
in both figures, r are rather small at 62 = 0.64, indicating that occupation 
of these intruder orbitals are not important for this deformation. However, at 
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Fig. 3. Percentage of high-j intruder occupation for different spin states with de- 
formations around the HD minimum in 110 Pd. 

£2 = 0.74 where the HD minimum is about to develop, a significant increase 
in r is observed. In addition, r increase rapidly with spin, suggesting that 
the highest j intruder orbitals play a more important role in high spin states. 
Finally at e<i = 0.84 beyond the HD minimum, r reaches the maxima and stay 
constant with spin (about 14.5% for the ii 3 / 2 protons, and 16.5% for the ji 5 / 2 
neutrons). Thus roughly, the highest j intruder contributes about 1/7 to the 
states under discussion. This is a remarkable value, considering the fact that 
this amount comes from only a few intruder orbitals. On the other hand, the 
major part of the contribution is indeed from many particles of the rest orbits 
although each individual may only contribute a little. 

Pairing is one of the most important correlations in nuclear states. In many 
high-spin calculations, however, pairing is often neglected for simplicity. By 
using the total wavefunction of Eq. (j2J), we have evaluated the pairing corre- 
lation in terms of pairing gaps 



Al 



(7) 



In Eq. (J7J), P T is the pair operator, G T M the pairing strengths in Eq. (J5J), and 
t = 7r or v. The results for the HD states, with the corresponding wavefunc- 
tions yielding a quadrupole moment (Q) = 3.9 eb, are plotted in Fig. 4. It can 
be seen that the pairing gaps decrease with spin as anticipated, which has been 
understood as a reduction in pairing correlation due to the rotational align- 
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Fig. 4. Calculated (proton and neutron) pairing gaps in the HD states in 110 Pd. 

ment. However, our many-body wavefunction does not give a sudden collapse 
in the pair field, or a phase transition from a superfluid to a normal state. 
This means that our results suggest that a considerable amount of pairing 
correlation remains even in the highest spin states in the HD well. 



4 Quantitative predictions 

We finally present quantitative predictions for yrast energies of the rotational 
bands (the lowest state at each angular momentum) in 110 Pd at each of the 
minima. In Fig. 5, 7-ray energies are calculated by the PSM with the basis 
deformation e 2 = 0.25 for the ND band, e 2 = 0.57 for the SD band, and 
e 2 = 0.79 for the HD band, respectively. The experimental ND yrast band in 
this nucleus is also shown for comparison. To support the prediction, the PSM 
results for the two experimentally-known bands in 108 Cd are also calculated 
by using the same calculation condition. As one can see, our results agree 
well with the known ND band in 110 Pd, as well as with the known ND and 
the highly deformed band in 108 Cd. The SD and HD bands in 110 Pd are our 
predictions, awaiting experimental confirmation. 

We have also carried out systematical calculations for the Zr, Mo, Ru, Pd, and 
Cd isotopes. Our results indicate that some isotopes of Zr, Mo, Ru, Pd, and 
Cd, with the neutron number around 64, show separate minima corresponding 
to SD and HD shapes. For all these nuclei, a SD minimum is found at e 2 ~ 
0.56 starting from I pa 24, and a HD one at e 2 ~ 0.80 from I pa 44. These 
predictions are obtained based on a theory that has quantitatively described 
the observed band in 108 Cd [91TTT] . which, according to our present results, lies 
just a step away from the mass region exhibiting hyper- deformed shapes. 
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Fig. 5. Calculated 7-ray energies (E^(I) = E(I) — E(I — 2)) for 110 Pd (upper panel) 
at different energy minima shown in Fig. 1. The theoretical results are compared 
with available ND data. Those for SD and HD states in 110 Pd are predictions. 
Results for 108 Cd (lower panel) are also given, and compared with known data. 



5 summary 



To summarize, using the projection method implemented in the Projected 
Shell Model, we have performed energy surface calculations for the A- 110 mass 
region. Our results for 110 Pd have indicated that three distinct energy minima, 
occurring at different angular-momentum regions, coexist in one nucleus. We 
have further investigated the source for the occurrence of these minima. The 
exact treatment of nuclear rotation through angular momentum projection 
can be important because the shapes in this study appear as sensitive func- 
tions of angular momentum (e.g. the rapid shape change in low -spin states at 
the ND minimum and the emergence of SD and HD minimum at sufficiently 
high-spins). Pairing correlation is sustained in our wavefunctions even at the 
highest spin of the HD states. Detailed analysis has shown that, while the 
nucleons occupying the high-j intruder orbitals have a larger amount of indi- 
vidual contributions to the total collectivity, the major part of collectivity is 
originated from the rest normal orbitals. 

Based on the present study, we have suggested that several nuclei in the A- 
110 mass region, with the neutron number around 64, would be the promising 
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candidates in which long-anticipated hyper-deformed rotational bands may be 
observed. Thus, the present work has shown a picture of coexistence of three 
distinct deformation minima, corresponding to normal, super-, and hyper- 
deformed shapes, in nuclei. We hope that these predictions can be experimen- 
tally confirmed in the near future. 

This work was supported partially by the Chinese Major State Basic Research 
Development Program through grant 2007CB815005, by U.S. NSF under con- 
tract PHY-0216783 (Y.S.), by U.S. DOE under contract DE-FG02-96ER40983 
(J.-y.Z.), and by NNSF of China under contract 10325521 (G.-L.L.). 
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